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Abstract 

The interaction of the electromagnetic field with a two dimensional plasma 
sheet intended to describe the pi-electrons of a carbon nano-tube or a Cgo 
molecule is investigated. By integrating out first the displacement field of the 
plasma or first the electromagnetic field different representations for quantities 
like the Casimir energy are derived which are shown to be consistent with one 
another. Starting from the covariant gauge for the electromagnetic field it is 
shown that the matching conditions to which the presence of the plasma sheet 
can be reduced are different from the commonly used ones. The difference in 
the treatments does not show up in the Casimir force between two parallel 
sheets, but it is present in the Casimir-Polder force between a charge or a 
neutral atom and a sheet. At once, since the plasma sheet is a regularization 
of the conductor boundary conditions, this sheds light on the difference in 
physics found earlier in the realization of conductor boundary conditions as 
'thin' or 'thick' boundary conditions in Phys.Rev.D70(2004)085010. 

1 Introduction 

The boundary conditions the electromagnetic field has to obey on the surface of a 
conductor, 

^11 = 5± = 0, (1) 

belong to the standard topics in classical electrodynamics. However, in [Ij it was 
shown that the conditions ([1]) do not fix the electromagnetic field completely, rather 
they allow for a freedom in the normal component of the electric field, on 
the surface of the conductor. In the standard treatment, which is mostly done in 
Coulomb gauge, this freedom is fixed by applying Gauss's law, div£^ = 0. In [1] 
it was shown that a consistent scheme can be constructed in which £"_]_ may well 
stay without any condition on the surface. The difference between both treatments 
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is that in [T] the surface was assumed to be mathematically thin, i.e., purely two- 
dimensional, whereas in the standard treatment a conducting bulk behind the surface 
is assumed, from which the condition on the normal component might be deduced. 
For that reason, in [1] the standard realization was called 'thick' conductor and the 
new realization in which no additional condition on £^ are assumed was called 'thin' 
conductor. The question to which extend the 'thin' conductor is physical, could not 
be answered in a satisfactory way in [Ij due to the lack of a physical model for an 
infinitely thin conducting sheet. However, it is desirable to have such a model since 
in it was shown that 'thick' and 'thin' conductors result in different physics. For 
instance, for 'thin' conductors there is an interaction of classical charges across the 
surface which is absent for 'thick' conductor. Furthermore, the Casimir-Polder force 
between an atom and a surface comes out by 13 % smaller for 'thin' as compared to 
'thick' conductor. It is interesting to note that the Casimir force between two flat 
surfaces is the same for both. It should be noticed that the difference between both 
approaches is not in the different gauges chosen but in the boundary conditions. In 
it was shown that in the approach with 'thin' boundary conditions and covariant 
gauge the standard result can be reproduced if imposing an additional boundary 
condition on the normal component of the electric field which does not follow from 

In the present paper we answer the question about a physical realization of 
the 'thin' conductor. As such we consider a two dimensional plasma sheet and 
show that it delivers in the limit of its plasma frequency to infinity just the 'thin' 
conductor and, for instance, the reduced Casimir-Polder force. The two-dimensional 
plasma sheet plays a twofold role. On the one hand side it is a regularization of a 
ideally conducting surface (its plasma frequency VL plays the role of the regularization 
parameter) which is recovered in the limit Vt ^ oo. On the other hand side it is a 
model for the vr-electrons of a graphene sheet, a Cgo molecule or a carbon nano-tube. 
The interaction of an atom with such a object is not only a measurable quantity 
but it is at present of obvious practical interest. 

In the next section we start from the general formulation of a system consisting 
of the electromagnetic field and an electric fluid (plasma) confined to some two- 
dimensional surface S. The fluid is described by a displacement vector ^ and we use 
the formulation in terms of a functional integral representation which is, as usual, 
not the only possible, but the most convenient formulation. We work in covariant 
gauge which in general terms looks more natural even in connection with the non- 
relativistic dynamics of the fluid. In order to make the treatment more convincing 
we consider two different ways by first integrating out either the displacement field 
^ or the electromagnetic field. In the first case we obtain the photon propagator in 
the presence of the plasma sheet and in the second one we get the propagator of 
the displacement field. In both cases we consider the excitation spectrum, i.e., the 
plasmons. To simplify calculations we first consider a scalar model before entering 
the discussion of the photon polarizations. As applications we (re-) calculate the 
Casimir force between two plasma sheets and the Casimir-Polder force between an 
atom and one sheet. In the appendix we consider the excitations for a spherical 
plasma shell. 

Throughout the paper we use units with h = c= 1. 
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2 The plasma shell model and its interaction with 
the electromagnetic field 

The plasma shell model describes an electrically charged fluid confined to a two- 
dimensional surface S" by a displacement vector ^(x) (x G S) which is tangential to 
S. A immobile, overall electrically neutralizing background is assumed. This model 
was probably first considered in [2] to describe a layered electron gas and, later, in 
[3] to describe the plasmon oscillations in Ceo- In fact, this simple model describes 
quite well a number of properties of the vr-electrons in Ceo, carbon nano-tubes and 
objects alike With respect to the Casimir effect it was considered in a series of 
papers, [5] and successors. In [6] the interaction of two such sheets with emphasis 
on the role of the surface plasmons and in [7] the interaction of one such sheet with 
a dielectric half space and with a neutral atom were calculated. As for the past 
two papers it should be mentioned that there the standard realization of boundary 
conditions was taken as granted and that both papers started from the boundary 
conditions on the TE and TM modes and standard polarizations. As it will be shown 
below the results of [6] remain valid whereas that of [7] need to be reconsidered. At 
least the result obtained there for the interaction of an atom with a graphene, i.e., 
with a plasma sheet, is different from what we will find below in section 5. 
The action of the considered model reads 

+ j d^x A^{x){j^{x) + J^{x)), (2) 

where F^y = d^A^ — d^A^ is the usual field strength tensor (we use the notation 
= d/dx^) and ^ is the displacement field of the fiuid having a mass density m. 
We use X, y, ... to denote coordinates in the bulk and z', . . . for coordinates on the 
surface S and we adopt the usual summation convention that for repeated indices 
the sum over their range is assumed. So the first term in S, ([2]), is the usual action 
of the electromagnetic field in covariant gauge with gauge fixing parameter a, the 
second is the kinetic energy of the fiuid (the dot denotes the time derivative and 
we assume non-relativistic dynamics) and the third describes the interaction of the 
electromagnetic field with the fiuid current j^{x) and some external current J^{x). 
The fiuid current is 



jo(a;) = e dz6{x-f{z))V^{z), 
Js 

j(x) = e I dz5ix-f{z))iiz), (3) 



s 



where e is the charge density. The surface is described by S* = {x | x = f{z)} and 
2; is a parameterization of S. In the simplest case of a plane perpendicular to the 
xs-axis at X3 = 0, a parameterization is z = {zo,zi,Z2} and f{z) = {zq, zi, Z2,0} 
defines the surface S. In that case the current is 



J,{x) = e6{x,){ >;/ (4) 
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and the displacement vector is ^{x) = { ^ {x) j ^'^ ~ 0)1)2). Current conserva- 
tion holds, d^i^ = 0. 

For the following it is convenient to represent the interaction in terms of an 
integral kernel, H^-i^b{x, z), 

d^x A^{x)j^{x) = j d^x J d^z A^{x)H^.i^i,{x,z)ii^b{z). (5) 

Obviously, this is possible for any surface 5* and it allows, in this way, to consider 
the most general case. However, here we restrict ourselves to the special case of flat 
surfaces, located at 0:3 = {h numbers the surfaces if more than one), 

Hi,-i,b{x, z) = eS{xa - Za)S{x3 - Of,) {Sf,odi - 5^ido) (6) 

(/i = 0, ...,3, a = 0,1,2, i = 1,2, b = 1,2). In the case of two surfaces the 
parameterization must include the number of the surface too. For a single surface 
the index b can be simply dropped. 

Now let us consider the generating functional of the Greens functions. 



Z{J) = j DAD^ e'^ (7) 
with the action S rewritten in the form 

S = ^Jd^xA^{x)K^,A,{x) + ^J^dh^,,b{z)K^^,,,{z) 

+ d*x d^z A^{x)Hij^.i^b{x,z)^i^b{z) + / (i^x y4^(x) J^(x), (8) 



s 



where 



1 



a 



= 9^..^' -ll--\d^d, (9) 



is the kernel of the free action of the electromagnetic field and 

= -md^ (10) 

is for the displacement field. 

The generating functional is in the known way by means of 

F=^lnZ(0) (11) 

connected with the free energy F (T is the total time and it drops out in a few steps). 
The distance dependent part of the free energy is then the vacuum or Casimir energy 
of the considered system. 

At this place a remark on the role of the displacement field ^ must be added. In 
letting its mass m go to zero the term with the kinetic energy disappears and the 
action S, Eq.([8]), becomes linear in ^. Then the functional integral over ^ can be car- 
ried out and it delivers a functional delta function, Yl ^(/ dx A^{x)Hfj,-i^b{x, z)). 



As it will be apparent below, this is nothing else that the treatment of conductor 
boundary conditions by functional delta functions introduced in [8] . In this sense 
the kinetic energy of the displacement field provides a regularization of these delta 
functions and of the conductor boundary conditions. 

In the functional integral ([7]), the action S, Eq.([H]), is quadratic in the fields, 
hence the integral is Gaussian and can be carried out. For the sake of clarity we 
adopt a symbolic notation dropping all indices and arguments. The the action is 

S = ^AKA + ^^K°^ + AH^ + AJ. (12) 

Now we first integrate out the displacement field ^. We write S as complete square 
in ^, 

S = ^A'KA + ^{^ + OK''{^ + + AJ, (13) 

with 

^={K'')-^H^A (14) 

and 

^K = K-H{K'')-^H^. (15) 

Here, is the transposed in the sense that AHE, = C,H^A must hold. The inverse 
[K'^)^^ of like similar inverse below must be taken in the space in which the 
operator acts, the surface S in this case. 

In a second step we complete the square for the field A in order to integrate out 
the electromagnetic field and represent the action in the form 

5 = 1 (A + i) (A + i) + 1 (^ + i) {i + A'K)-'J (16) 

with A = {^K)~^J. Now the Gaussian integrations can be carried out, first that 
over ^, subsequently that over A and we arrive at 

Z={det'Ky''^ {detK')-"^ e'^^^"^)-'^. (17) 

In these formulas, is the kernel of the action of the electromagnetic field after 
integrating out the displacement field, i.e., taking into account the interaction with 
the plasma sheet. In general, we denote by an index 5* in front a kernel that includes 
the interaction and by the index we mark a kernel belonging to the displacement 
field. 

From the action f|T3|) . the equations of motion for the field A follow, 

^KA = J, (18) 

where the interaction with the fluid is in the second term in the r.h.s. of ^K, Eq. ffTSl) . 
For vanishing source J = 0, these are the equations determining the modes of the 
electromagnetic field in the presence of the fluid. Below, in an example, we will 
see that provides the usual free space Maxwell equations supplemented by the 
matching conditions on the surface S. 
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In this formulation, the Casimir energy is given according to flTT]) by 

F = — trln^ir+ — trlnK°, (19) 

where the traces are taken in the corresponding spaces. For a setup with two sheets 
the distance dependent part is contained in the first term, i.e. in tr In ^K, whereas 
the second term, trlnK", adds only a (infinite) constant. 

Now we consider the case of first integrating out the electromagnetic field A. For 
this, we rewrite the action (fT2!) in the form 

S = \{A + A) K {A + A)+^-iK'' i-]^AK A, (20) 

where now 

A = K-^{Hi + J). (21) 

Here is the inversion of the kernel of the free action in the empty space, i.e., 
the free space photon propagator D = K~^, in detailed writing 

(+iO for the causal propagator). After that we complete the square for ^. By means 
of 

- AK A = -iH^K-^Hi + iH^K-^J + -JK'^J (23) 

we get 

S=^-{A + A)K{A + A) + ]^i'^K' i-iH'^K-'j-^-J K-^ J (24) 

with 

^K^ = K^ - H^K~^H. (25) 

Obviously, ^K'^ is the kernel of the action of the displacement field after integrating 
out the electromagnetic field and the equations of motion for the fluid read 

^K^ ^ = (26) 

(below we will investigate this equation in detail). The first term in '^K^, Eq. fl25l) . 
results from the kinetic energy and the second from the electromagnetic interaction 
of the fluid with itself, for instance, due to the Coulomb interaction. 
Completing now in Eq. (!2^ the quadratic form for ^ we arrive at 

S=^iA + A) K iA + A) + ^i^ + i) ^K' {^ + i)-}-jSDj (27) 

with i = {^K^)-^H^K-^J and 

= D - DH(^K°)-^H^D, (28) 
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which is the photon propagator in the presence of the fluid. Now the Gaussian 
integrations in ([7]) can be carried out using the action (1271) . first that over the A- 
field, and then that over the ^-field, 

Z={detK)-"^ (det V)"'^' e-^''^'. (29) 
The corresponding Casimir energy is 

F = — trlnK + — trln'^ir". (30) 
2T 2T ^ ' 

Thanks to the derivation, this expression must be a representation of the same 
Casimir energy as in (1T9|) . The distance dependent part in (130|) is in the second 
term, whereas the first is distance independent. 

As seen, (1191) and (l3Up are two different representations of the same quantity. 
This follows from the derivation. For instance, if taking a finite dimensional example 
where the kernels are represented by finite dimensional matrices acting in different 
spaces, this equivalence can be shown quite easily. In field theory where these 
kernels are infinite dimensional matrices, it is not so easy mainly due to the inherent 
divergences. However, there is no doubt that correct definitions can be given and 
that the equivalence can be shown. We are not going to do that and restrict ourselves 
to some examples. 

First of all, let us compare the source terms in (1T71) and (l29l) . Since they must 
be equal, 

^K~^ = ^D, (31) 

is the inversion of the kernel of the electromagnetic field after integrating out 
the displacement field. In this form describes the modification of the photon 
propagator which comes in from the surface S. In fact, an expression of this type 
was derived in [S] in the simpler case of conductor boundary conditions and in 
[H Uni for delta-function potentials (later it was re-obtained by other authors and, 
for instance in [Tl] it was discussed in connection with the Krein formula). 

Now we consider in the l.h.s. of Eq. (l3T]) in the form given by Eq. dlSj) . Up 
to the Lorentz indices it consists of the wave operator and a second term which 
due to the delta functions in j^, Eq.(I3]) or (jlj), has support only on the surface S. 
It can be considered as a potential with localized support and in the simplest case 
we get a wave equation with a delta function potential like the one considered in 
ini HOj . Since such potentials are equivalent to known matching conditions the 
electromagnetic field has to obey on the surface the kernel is equivalent to 
the free wave equation with the corresponding boundary conditions. In this way 
Eq.(pTl) can be verified. 

Now we turn to the discussions of the expressions (1T9|) and (130|) for the Casimir 
energy and consider, say for two plasma sheets, the distance dependent part. In 
( 1T9|) it results from and in (130|) from ^K^. We note that these are different 
operators, acting in different spaces, in the bulk and ^K*^ on the surface S. Also 
their physical interpretation is different, '%'describes the photon fluctuation in the 
presence of the surface 5* and ^K^ describes the fluctuations of the displacement field 
^ including its electromagnetic self-interaction. The contribution from both to the 
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Casimir energy must be the same. Moreover, this contribution enters only one time, 
either through in representation (fT9!) or through ^K^ in representation (130|) . but 
not twice. This is quite counterintuitive since one is tempted to argue that both 
kinds of fluctuations should contribute to the vacuum energy. But since that is not 
the case one can conclude that in a scheme of canonical quantization it would not 
be possible to introduce independent creation and annihilation operators to both 
kinds of excitations. 



3 A scalar field interacting with plane plasma sheets 

In order to have an explicit realization of the general formulas developed in the pre- 
ceding section we consider the simplest non-trivial example, a scalar field interacting 
with one or two parallel plasma sheets. We start from the action ([8]) and replace 
the electromagnetic field A by a real scalar field 0, 

+ f d'x <j){x) Uoix) + J{x)) . (32) 



The kernel of the free action of (p is now simply the wave operator, 

K = -d'^ = -d'^ + A , (33) 

K'^ = —mdQ is the analogous quantity for the displacement field and the current is 
given by 

j{x) = e ^ S{x3 - ab)Vi^ifi{x). (34) 

6=1,2 

In the case of two plasma sheets these are located at x^ = at {b = 1,2). In this 
model, e can be viewed simply as a coupling constant without having the meaning 
of an electric charge and we define 

Q=- (35) 
m 

as plasma frequency. Strictly speaking, e and m are the corresponding densities per 
unit area of the sheets, however, since only the ratio fl35|) enters the final formulas 
we do not need to introduce separate notations for the densities. The displacement 
of the fluid is described by a two dimensional vector ^i^b{xa) in the plane of the 
sheet. Its components are labeled by i = 1, 2, Xa with a = 0, 1, 2 are the coordinates 
in the sheet and 6 = 1,2 labels the sheets. The interaction between the field and 
the displacement field ^ can be rewritten in the form of Eq.([6]), 



d^x (f){x)j{x) = I d'^x I d^z (j){x)Hi^i){x, z)^i^b{z) (36) 
with 

Hi^b{x, z) = e6{xa - Za)S{x3 - ab)Vi . (37) 
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The transposition of H, 

Hjh{x, z) = -e5{xa - Za)5{x3 - ab)Vi (3^ 

follows from 



d'^x J d-^z (f){x)Hi^b{x, z)^i^b{z) = j d^x J d^z ^i^b{z)Hlf,{x, z)(j){x) (39) 

with integrating by parts. 

For the simple geometry of two parallel sheets all these quantities can be sim- 
plified by Fourier transform in the x^-directions {a = 0,1,2). We define for the 
fields 

/d^k 
J^3 e^'--^ M^3) (40) 

and for the kernels 

K{x,y) = e^'='^(--^") K,{xs,y3) (41) 

and similar for ^ and The transformed quantities are marked by an index k. 
Since the fields are real, their Fourier transforms obey (pki^s)* = (p-kix^). The wave 
operator is connected with the corresponding kernel by K{x,y) = 6{x — y)K{x) and 
after Fourier transform it becomes 

Kkixs) = T' + dl (42) 

with 

r = Vfc2 + tO = ^kl-kl-kl + iQ , (43) 

where we introduced the infinitesimal imaginary part for later use. The correspond- 
ing quantity for the displacement field is simply 

Kl = mkl (44) 

and for the interaction kernel we note 

Hki^b = 'i'<^kiSix3- ttb). (45) 

Another simplification occurs if introducing polarizations for the displacement field. 
Obviously there are two polarizations. 



We use the notations TE and TM because the in the electromagnetic case these 
modes couple just to the corresponding polarizations. From the explicit form (H5|l 
it is clear that because of 

Hki,b^i,b = iekii6{x3 - ab)^l'^ (47) 
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only the TE-mode of the displacement field couples to the scalar field. The effective 
coupling constant is ek\\. 

Now we integrate out the displacement field ^. The first quantity we have to 
consider is the inversion of entering Eq. (fT5|) . Using the Fourier transform this is 
simply 

with 

= . ■ (49) 

^ mkl + iQ ^ ' 

Here and in the following we always chose the causal propagator. Applying H from 
the right and from the left and inserting into (|T5|) and using ^K^^x^, 2/3)^(3^3 — ^3) = 
^Kkix'i) we get 

'K,{x,) = T' + - ^(^3 - a,), (50) 

where k\\ = a/^J + /c| is the momentum parallel to the sheets. 

Now we consider the equations of motion for the field which are analogous to 
Eq. ffTSl) . Without source these read 



and, rewritten, 



''Kkixs) M^s) = (51) 



+ fc2^ E '^(^s - at) ] M^s) = rVfc(x3). (52) 



6=1,2 



Obviously, this is nothing else than a Schrodinger equation with delta function 
potential, a favorite textbook example. In the given case the potential is repulsive 
(f2 > follows from fl35|) independently on whether e has an interpretation as electric 
charge or not) and the field has a continuous spectrum. 

It is well known that the delta-function potential can be reformulated in terms of 
matching conditions on the surface. The field must be continuous and its derivative 
has a jump, 

Qkf, 

discont (f)' = -jY- (p ■ (53) 

In quantum field theory such potentials were considered several times. We refer 
here to [HI [ID] , where delta potentials on two parallel planes were considered for 
the scalar field and also for a spinor field. There the propagators were written down 
and, using the energy-momentum tensor, the Casimir force between two planes was 
calculated. The latter is equivalent to calculate the Casimir force using Eq. (fT9|) so 
that we do not need to repeat that calculation here. 

Now we consider the second approach where the scalar field is integrated out 
first. The Fourier transform of the free space propagator follows by performing the 
integration over k^ in fl22|) (dropping the bracket containing the vector structure), 

gir 1x3-2/3 1 

Dk{xs-y3) = ^^ (54) 
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and from fl2^ we get 

{'K',U,,, = - e^r|a.-vl (55) 

where we added the necessary indices. This kernel is defined on the surface S, hence 
after Fourier transform it is purely algebraic, i.e., it does not contain any derivative 
or dependence on the coordinate X3. The corresponding equation of motion for the 
displacement field is given by Eq. (12Bi) . It describes the excitations of the plasma 
including its selfinteraction caused by the scalar field. To consider these excitations 
it is sufficient to consider one sheet. In that case and after division by m the 
equations of motion read 

klS.^ - = . (56) 

These are diagonalized by the polarizations fj46l) and the eigenvalue equations be- 
come 

^o-^lf = 0, (TE) 

kl = 0. (TM) (57) 

The meaning of the second equation is simply that this mode does not couple to the 
scalar field and it does not have excitations. The equation for the first mode has no 
solution for > (taking into account that F has a positive imaginary part). 

Now we consider the construction of the propagator for the 0-field which is 
given by Eq. fl25]) in general terms. In fact this repeats the corresponding calculation 
in P, [To] . Therefore we restrict ourselves here to the simplest case of having only 
one plane. In that case, after Fourier transform it becomes 

^Dfc(x3,i/3) = — ieki (^K^)r^ iek. — — (58) 

with the inversion 



U ZlL 

which is done using the polarizations from (j46l) . The second polarization does not 
contribute and we get finally 

„ir\x3-y3\ gjr(|a;3| + |;/3|) 

^D,{x^,y^) = ^^ — r(fc), (60) 



where 



= (61) 

^ ^ 1 2jr fco 



n k 



has the meaning of being the refiection coefficient in the corresponding scattering 
problem. This propagator (!60!) obeys the matching conditions (!53|) and up to nota- 
tions it coincides with the scalar propagator derived in ^ (TU] . 

In this way we demonstrated how the general formulas of section 2 can be realized 
in the simplest case of a scalar field interacting with plane plasma sheets. Especially 
the equivalence of different representations was demonstrated. 
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4 The electromagnetic field interacting with plasma 
sheets 



In this section we consider the electromagnetic field interacting with one or two fiat 
plasma sheets. On the one hand side this is a straight forward generalization of the 
scalar case considered in detail in the preceding section, but on the other hand side 
from the electromagnetic case the most interesting conclusions follow. 

We start with introducing polarizations for the electromagnetic field. We do not 
chose the standard polarizations but those first introduced in [12] (and later in [13] 
for a sphere and a cylinder). After Fourier transform ( l40l) the polarizations vectors 
read 



h 

\ / 






V 1 / 



K 



( \ 

-A:i 

V / 



, El 



( k\\ 

koki 
kok2 




(62) 



These form a basis, 



g,, = gst K (63) 

with Qst = diag(l, —1, — , 1 — 1). Note that these vectors do not contain k^ and that 
in the case of conductor boundary conditions only the polarizations with s = 1, 2 
are affected by the boundary. In the approach of the 'thin' conductor, those with 
s = 0, 3 were kept free of conditions. We will see that this takes place also in the 
case of the plasma sheet. 

With fl62|) and the Fourier transform (140|) . the expansion of the electromagnetic 
potentials reads 

s=0 



(64) 



where we showed the sum over s explicitly. We note that without the plasma sheet 
the equations of motion for the amplitudes are (F^ + ^K^^s) = with F defined 
in fH3|) and that the free part of the action becomes diagonalized when inserting f lM|) 
into (IHl). We start from inserting ( IMl) into the interaction with the plasma sheet, 
i.e., into ([5]). After Fourier transform we note (for one plasma sheet) 



{Hk)^,,i = ie6{x3) {6^oki - d^iko) 
and the together with the polarization vectors we get 

0, (s = 0,3) 



(65) 



hi 



k ) 



-k2 

ki 

ki 
k2 
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2) 



(66) 
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In this way, only the polarizations fl62l) with s = 1, 2 couple to the displacement 
field of the plasma sheet. This is the most important observation since it allows 
to diagonalize the action completely and to establish the correspondence with the 
'thin' conductor approach in yj . Using the mode expansion (H6l) of the displacement 
field we get for the interaction (jS]) 

J d'x A,{x)Ux) = / ^ {ko 4™ AliOy - T 4™ ^^(0)*) • (67) 

Here Al{0)* (s = 1, 2) are the complex conjugated amplitudes taken at the position 
of the sheet, i.e., at = 0. In this way the action is diagonalized and the interaction 
of the electromagnetic field with the plasma sheet is reduced to two scalar problems. 
More exactly, the complete action becomes a sum of thee parts. The first contains 
the polarizations with s = 0,3. This part does not depend on the plasma sheet. 
It contains the gauge fixing parameter a. The next two parts, with s = 1 and 
s = 2, are just the two scalar problems where the plasma sheet enters and these 
are independent on the gauge fixing parameter a. For them we can use the results 
of the preceding section with the substitution of the effective coupling constants 
introduced in Eq. fHTl) according to 

efc|| efco, (TE) 

efcii er. (TM) (68) 

These substitutions follow from comparing Eq. (HSj) with (16 7p . In order to use them 
for the photon propagator we represent it in the form 

'D,,{x,y) = Dl:='''\x,y) + ^ ^^^(x,!/). (69) 

s=l,2 

Here the first part carries the gauge dependence and it is independent on the plasma 
sheet. The remaining two parts read in Fourier representation 

(s = 1,2, no sum over s here) with 

, , pjri^s-ysi p«r(|x3|+|j/3|) 

= — 

= D['>(x3.m)-Di"{x3,y3), (71) 

where in the last line the separation into free part and boundary dependent adden- 
dum was introduced which we will use in the next section. The refiection coefficients 
rs{k) contain the different couplings (1681) relative to the scalar case which can be 
accounted for by the formal substitutions of the plasma frequencies 

{s = l) 

= { ' (72) 
n^. is = 2) 
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With these relations, the reflection coefficients are 



ri(fc) 



1 



2ir ' 

n 

1 



(73) 




From a comparison with [14J it follows that ri{k) is just the reflection coefficient for 
the TE-mode which is completely obvious and that r2{k) is the reflection coefficient 
for the TM-mode which becomes obvious after the remark the by virtue of the 
equations of motion F is the same as a momentum perpendicular to the sheet if one 
introduces such. 

The difference between the treatments of the scalar problems corresponding to 
the TM-mode here and in the standard approach is the following. The match- 
ing condition for the TM-mode in the approach taken here implies that the mode 
function is continuous across the sheet and that its derivative has a jump given by 
Eq. (153!) with Q substituted by Qg from (172!) and in this way it corresponds to a delta 
function potential. In the standard approach one introduces another mode function 
whose derivative is continuous and whose value across the sheet has a jump. That 
corresponds to a potential given by the derivative of a delta function. As for the 
resulting scalar problems, these are obviously the same as far as they depend only on 
the reflection coefficient r2{k), Eq. fl73|) . From here as a consequence it follows that 
the Casimir force between two sheets is the same in both approaches. Remember 
that it can be calculated by a generalized Lifshitz formula where only the reflection 
coefficients enter. This is also the reason why the results of [B] remain valid. 

However, this does not imply that the Casimir force between a plasma sheet and 
a dielectric half space must be the same in both approaches since for the dielectric 
half space one would use the standard polarization and for the plasma sheet the new 
ones. 

In the standard approach it is known that for the TM mode there is a surface 
plasmon and that in the TE mode there is none. The same hold in the present 
approach too. For the TE mode the matching condition follows from fl53|) with the 
upper line of fl72l) and it reads simply 



where (p stands for AKx^) introduced in Eq. flMl) . This corresponds to a repulsive 
delta potential. In a similar way, for the TM mode we get with the lower line in 



where stands for A^xs) introduced in Eq. (!Ml) . As can be seen, for /cq < fc^ the 
delta potential becomes attractive. In that case the one-dimensional problem has 
one bound state and this bound state is a wave propagating along the surface in 
the three dimensional case. These surface plasmons were considered for instance in 
[T4] . It should be mentioned that these are to a large extend analogs of the surface 
plasmons known to travel on a flat surface of a dielectric body if for the permittivity 



discont (f)' = Q (j) , 



(74) 



discont (f)' = Q—^ (p , 



F2 



(75) 
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the plasma model is used. The same statements can be obtained from the behavior 
of the reflection coeflicients (173|l on the upper imaginary momentum axis. 

There is another way to consider these waves. Namely, let us consider the ex- 
citations of the plasma, i.e., the excitations of the displacement fleld. These are 
subject to the general equation (12^ . Knowing that we have two scalar problems 
for the amplitudes ^™ and ^™ of the displacement fleld, which were introduced in 
Eq. fH6l) . with plasma frequencies given by Eq. (l72l) . we get from the corresponding 
scalar equation, namely from the upper line in fl57j) . the equations 

0, (TE) 

0. (TM) (76) 

For the TE-case we note that for > because of the positive imaginary part of 
r we do not have a solution besides the trivial one ko = 0. In opposite, for the 
TM-case we have a solution for fcg < A;^, 

kl = ^^^' + ^QH-Y- ^^^^ 

We note that this spectrum is just the same as for the surface plasmon following 
from Eq. flTS]) mentioned above. There is reason to think that this coincidence is not 
accidental. Let is return to section 2. The equations of motion for the displace- 
ment fleld are given by Eq. (l26|) . After diagonalization, the zeros of ^K^ deflne the 
spectrum. On the other hand side, the spectrum of the electromagnetic fleld can 
be found from the poles of its propagator which in the presence of the plasma 
sheet is given by Eq. (l28!l . There in the right hand side we have ^K'^ in the denom- 
inator so that it zeros deflne poles of ^D. In this way one can show in the general 
case that the surface plasmons and the excitations of the plasma (which are usually 
called plasmons) have the same spectrum, or more exactly, that for each mode of 
the displacement excitations there is a corresponding mode in the electromagnetic 
spectrum. In the Appendix we illustrate this on the example of a spherical plasma 
sheet. 



° 2i 



5 Interaction of a charge with a plasma sheet 

In this section we consider the interaction of a charge and of a neutral atom with a 
plane plasma sheet. For the neutral atom this is the same setup as for the Casimir- 
Polder force in [15]. We use also the same method, the quantum mechanical per- 
turbation theory. Since the original paper, the Casimir-Polder force was considered 
repeatedly, see for instance [16] or [T7|. In all cases, for the photon propagator the 
same setup was taken which in [T] was called 'thick' boundary condition. In [T], the 
Casimir-Polder force was calculated for the 'thin' boundary conditions and found to 
be smaller. In this section we repeat that calculation for a plasma sheet in place of 
the conducting surface used in [T]. It will be shown that we reproduce the result for 
'thin' boundary conditions in the limit of inflnite plasma frequency fl. 
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The calculation of the Casimir-Polder force is quite tedious and involves a num- 
ber of subtle moments. All these are known in literature and we do not need to 
repeat them here. Instead, we follow very closely the derivation given in [1] and 
focus on the modifications which come in from the plasma sheet in place of the 
conducting surface. In the initial setup of the problem we have three systems, the 
electromagnetic field, the displacement field and the charge resp. the atom whose 
interaction with the plasma sheet is the subject of our interest. We eliminate the 
displacement field by, say integrating it out in a functional integral setup like in sec- 
tion 2. In this way we assume that there is no other interaction between the charge 
and the plasma sheet than through the electromagnetic field. After that we we are 
left with the photon propagator in the presence of the sheet, i.e., obeying match- 
ing conditions as discussed in the preceding section. This propagator is given by 
Eg. (128!) in general terms and by Eqs. (l69!) to (17T]) for a fiat sheet. With this photon 
propagator we repeat the calculation done in section 2 in [I]. We note for instance 
that we use here the same polarizations E'^ and that the whole difference is in the 
presence of the reflection coefficients ri 2{k), Eq. fl73p . which are now present in the 
photon propagator. For the conducting surface, of course, ri 2 = 1 holds which for 
the plasma sheet is recovered in the limit Q — > 00. 

The calculation starts from a quantum mechanical Hamilton operator, Eq.(26) 
inP, 

(p - eX) 

H = ^- ^ + eAo + V{x), (78) 

2m 

for a particle (electron) with momentum p in a potential (of the nucleus) V{x) and 
interacting with the electromagnetic potential An{x). We note that e and m in (!78|) 
do not need to be the same as for the plasma sheet whose properties enter only 
through its plasma frequency fl, (135|) . 

The first quantity to be calculated is the static interaction of the charge with the 
nucleus, or more exactly, the change of this interaction brought in by the plasma 
sheet. In [1] this was AV, Eq.(31), which could there be written down directly as 
the potential of the mirror charge of the nucleus. In general it is given as the zeroth 
component of the field generated by a static source. In terms of the propagator this 
is 

AV{x) = e'J Doo{x,y)6{y) , (79) 

where we included only the boundary dependent part DQo{x,y) of the propagator 
which in our case is given by Eq.( 17T]) . In momentum representation, the integration 
over yo makes the momentum ko vanish and with formulas f l69l) to ( 17Ti) we obtain 

f dk\i kf, ir(2a-x3) 

Here we have taken into account that following [1] the sheet is located at X3 = a 
and the atom in the origin. A contribution to the (0, 0) component in the boundary 
dependent part comes only from the polarization vector E'^ in ( !62l) and it delivers 
the factor A;^/r^. Further we have to take into account that for fco = we have 
F = ik\\ so that this factor simply becomes (-1). The refiection coefficient r2{k), 
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given by Eq. fl73p simplifies too, for = it becomes r2{kQ = 0) = 1. Tliat means 
we get tlie same result as for the ideal conductor. In fact, performing the integration 
in (IHOj) we just get 



AV{x) 



and further, as a first order perturbation of the static energjl^, 

4.^„„, = (n|AV'(x)|„) = |l(i- + ^ + ...) (82) 

follows, where Q is the quadruple moment. 

The next step is to calculate the corrections following from y4(x)^-term in the 
Hamilton operator. Its boundary dependent term 6i is given by Eq.(36) in [1], 

1 _ 

Si = 7 < n\Dii{x,x)\n > . (83) 

2m I 

We insert the propagator into this formula and pay attention to the contributions 
from the two polarizations. 



li^i / If ^ (^■'*' ^ " " I I " > ■ 

In this formula one needs to perform the Wick rotation, —>■ ik^, F — > Z7 = 
i^Jkl + A;^ and after expanding the matrix elements we get in generalization of Eq. 
(38) in py 

with 

1 1 

Since the integrations in this formula cannot be performed explicitly we represent 
6i in the form 

^1 = -f (fTEi^a) + \hM{na)) +o(^) , (87) 

where the functions 

k 



poo 

hE{x) = dk 
Jo 



1 + k/x 



roo 



fTuix) = 3x I dk e \ 1 — -^/^ arctan W — 



10 



X k 



"'^ Following the notations in [16] and [T], till the end of this section interaction energies are 
denoted by 6. 
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describe the reduction of the corresponding contribution due to the plasma sheet 
relative to the ideal conductor case. We note 



/te(oo) = /tm(oo) 



(89) 



and, for x ^ 0, 

fTE{x) ~ X + . . . , /TM(a;) ~ 3a; + . . . . (90) 

Next we consider the second order perturbation which in [T] is given by Eq.(41). 
Again, the change coming in from the plasma sheet is to insert the reflection co- 
efficients into the contributions from the corresponding polarizations. Proceeding 
directly to the expression (43) in |lj the second order perturbation reads now 



EE- 

n' s=l 



1 f d^ka -e^'^" 



(27r)3 -2iT 



r,{k) 



\< n\ Gs\ n' >\'^ 
-ko + En- En' {I - ie) 



(91) 



For the matrix elements the expressions (43) and (44) in [IJ remain valid. 

The calculation advances with considering separately the cases of a single charge 
and a neutral atom in front of the sheet. For a single charge one uses the so called 
'no-recoil' approximation which implies an expansion for small energy differences 
En — En' in the denominator. 



— kn + En — En' 



1 



En — En' {En — En')" 



^2 



h.3 



+ 



(92) 



Using the expansion of the matrix elements given in [T] and the sum rules listed in 
[1] in the appendix we come to 



rno— recoil 



d^k^ -e2*r» -1 



(27r)3 -2iT ko 

,2 
72 



e^r2{k) -1 + 



e 

+ — 



kl k 
n(A:)+r2(A;)l-| + 



r2 ' klv^ 



+ 



(93) 



Now the integration can be carried out as described by Eqs.(49) in [T] as pole 
contribution delivering 



rno— recoil 
^2 



dk\\ 



g-2A:||a 



2k\\ 



— e 



e2 

+ — 



1 + ^ 



+ 



— e 



57ra 



-ak\\ 
1 



, , 2k\\ 

-'"^" + 2 + 7r 

2\ 



1 

2 



k 



ipi) 



(94) 
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with 

= / dke~'' ~] , +^ + ^ + - 
io \l + k/x 2 2 X 

hsix) = 1 + - . (95) 

X 

This is the interaction between a single charge and the plasma shell. The first 
term in flMl) is the electrostatic contribution which does not depend on the plasma 
frequency fl and it coincides with the ideal conductor limit. The second contribution 
depends on Q. For f2 — > cxo it turns into the ideal conductor case calculated in [Ij for 
the 'thin' boundary conditions. The weakening of this interaction energy relative to 
the conductor case is described by the functions h\\{x) and h^lx) which both obey 
h\\{oo) = h^loo) = 1. In the opposite limit, i.e., for small Q, these functions diverge 
which means that this limit is in conflict with the approximations made in deriving 

We note that the interaction between a single charge and a plasma sheet was 
calculated in [18], section 4, within the standard approach using Coulomb gauge. The 
result is different from (|94l) in the part proportional to (p^) (the contribution pro- 
portional to is the same). This result was to be expected since in [IB] the ideal 
conductor limit is what in [T] is called the 'thick' conductor. It must be mentioned 
that, as discussed in [1], the 'thin' conductor boundary conditions were derived for 
an infinitesimal thin conducting sheet. The question whether these apply to a con- 
ducting surface with a conducting bulk behind is still opened. It can be expected 
that in that case rather the 'thick' conductor boundary conditions apply. How- 
ever, this question needs further investigation. In opposite, the interaction with the 
plasma sheet, as follows from the above calculations, is uniquely described by flMl) . 

Now we turn to a neutral atom in front of the plasma sheet. The starting point 
is Eq. fl9T]) . This integral is split into two parts according to Eq.(51) in [1]. To the 
first part, which is the pole contribution in = 0, only the intermediate state 
with n = n' contributes. Since the momentum operator p in the matrix elements 
(see Eqs.(43) and (44) in [T]) has no diagonal contribution we have a non- vanishing 
contribution from the second polarization only. But the reflection coefficient r2{k), 
(173|l . is for ko = equal to unity. Hence this contribution is the same as for ideal 
conductor, ie.e, it is the same as calculated in p], Eq.(55), and it cancels against 
other static contributions. The second part in 0911) results from the VP-integral (see 
Eq. (51) in [T]) and after performing the Wick rotation it becomes 



rVP 



^ ' n' s=l,2 



where are given by Eq. (!86l) . In the next step we expand the denominator for 
small ^4 and use again the sum rules listed in the appendix in [1]. After that 
can be represented in the form 




hik)kl + r.ik)^ 1 - r,{k)kl I , (97) 
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where is the static polarizabihty of the atom in the i-th direction. Finally we use 
spherical coordinates in the ^-integrations with e = cos^, 



-e (he-' + f, (e' + (1 - ,r-)) - h(l - «")y } ■ (9i 



Here the first contribution cancels against 6i, flHTl) . and the rest is what contributes 
to the Casimir-Polder force. It can be written in the form 

5cp = I (^?TE(^^a) + ^gTui^a)^ + gsi^a) a,^ , (99) 

where the functions 

gTE{x) = 7: / afc 



6 Jo 1 + k/x 



5 r ,3 -k . + fl - e^^^ 



fi'TM(a^) = — / dk k e de 



22 Jo Jo 1 + e^k/x ' 

5 /""^ „ , s t /^2x 2x(x + k) v^(2a;2 + 2xk + 
— / dkk'^ e " \ — — — - + - — ^ —7;^ arctan ■ 



1 



5f3(x) = - dk k'^ e " de 



1 1 .2 



4 io Jo 1 + e^^/a; 

1 P „ , 3 . / a; Jx(x + k) [k\ , , 

describe the decrease of the corresponding contributions relative to the ideal con- 
ductor case considered in [IJ. All these functions are normalized to unity for x ^ 00. 
Also, they vanish for i7 — > so that within the approximations made the interaction 
with the plasma sheet disappears in this limit. 

The interaction of an atom with a plasma sheet was considered in [7J , Eq.(25), 
following the standard treatment and, consequently, it comes out to be different 
from ([99]). 



6 Conclusions 

In the foregoing sections we considered a plasma sheet interacting with the electro- 
magnetic field. We used a representation in terms of a functional integral. From 
that it was seen that the fiuid acts as a regularization for the conductor boundary 
conditions which are obtained in the limit of ^2 ^ 00, where Q is the plasma fre- 
quency of the fiuid. The consideration of this plasma shell model was motivated by 
the wish to have a physical model for the realization of the 'thin' boundary con- 
ditions introduced in [1] for a infinitely thin conducting sheet. Indeed, the plasma 
shell model may serve as such since it is a meaningful model for the vr-electrons of 
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a graphene sheet. For the electromagnetic field we used the same polarizations E^, 
( !62|) . as for the 'thin' boundary conditions in [1]. From this and from the refiection 
coefficients ( !73|) . it is clear that in the limit — >■ cx) we get just the 'thin' boundary 
conditions. It should be underlined that this is insofar a nontrivial result as this 
derivation is straightforward and it does not leave room for any freedom or arbi- 
trariness. This follows from the procedure applied consisting in first integrating out 
the displacement field which is done without any assumption on the gauge fixing 
or the polarizations of the electromagnetic field. After that we have an action for 
the electromagnetic field with the kernel ^K, Eq. ffTSj) and the corresponding equa- 
tions of motion, ( |T8l) . This kernel is uniquely defined and cannot be altered. The 
'remaining task', then, is the solution of the equations ( |T8l) . We use polarizations 
which diagonalize this kernel and come to two separated scalar problems. The 
one (TE-mode) is identical to the corresponding one in the standard approach, the 
other one (TM-mode) is different in that it has a different polarization vector (the 
corresponding reflection coefficients are the same). 

We would like to note that the attempt to solve the mentioned equations using 
the standard polarizations would hit the problem that these do not diagonalize 
as can be seen easily. Also, we do not see how the standard transition to Coulomb 
gauge can be done here since the equation of motion for the component Ao of 
the vector potential is not instantaneous (in opposite to the case when there is 
no displacement fleld). At the moment we do not have a good explanation what 
might be the problem within the standard approach. We restrict ourselves here to 
the repetition of the remark done already in the Introduction that in [1] we have 
shown that the standard approach, i.e., the 'thick' boundary conditions, follow from 
the 'thin' boundary conditions by introducing by hand a restriction on the normal 
component of the electric fleld in addition to ([T]). 

We would like to summarize the most sensible conclusions as follows: 

1. For an ideal conducting body, i.e., for a surface with boundary conditions 
^11 = B± = and a conducting bulk behind (this is the standard treatment, 
we call it 'thick boundary conditions'), our results do not apply. 

2. For a infinitely thin surface, the boundary conditions S\\ = B± = can be 
realized in a different way which we called 'thin boundary conditions'. These 
describe different physics, for instance, a by ~13% reduced Casimir- Polder 
force. 

3. For a two dimensional plasma sheet we get a unique description which is 
different from the standard one and which delivers, for instance, a different 
Casimir-Polder force. In the limit of inflnity plasma frequency when the sheet 
becomes a inflnitely thin conducting surface, we obtain the thin boundary 
conditions. 

As for the physical consequences in relation to the standard approach we showed 
that the Casimir force between two plasma sheets is the same whereas the Casimir- 
Polder force between an atom (or a charge) and a plasma sheet is different in both 
approaches. We cannot by-pass the conclusion that our approach gives the correct 
result. 
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In the second section we considered the calculation of the functional integral (it 
is Gaussian) in two ways, by first integrating out either the displacement field or 
the electromagnetic field. Different representations follow whose comparison allows 
for some insights. For instance, for the Casimir energy we obtained two different 
representations. In (HM the distance dependent part of the Casimir energy can 
be attributed to the vacuum fluctuations of the electromagnetic field whereas in 
f l30|) to those of the displacement field. In fact, this is a statement on what are 
the independent degrees of freedom, namely either the first or the second ones but 
not both at the same time. Furthermore, since the same Casimir energy can be 
calculated either from the modes living in the bulk or from modes living on the 
surface one could think of some kind of holographic principle. 

Another observation concerns the excitations. It turns out that the surface 
plasmons of the electromagnetic field have the same spectrum as the excitations 
of the fluid (which are usually considered the plasmons). We confirm the known 
results that there is no plasmon for the TE mode and that there is one (for each 
plasma shell) for the TM mode which has the same dispersion relation as known 
from [13]. In the Appendix we considered a spherical plasma sheet confirmed the 
corresponding conclusions (there is no plasmon at all). 

We have seen that there exist two realizations of the boundary conditions S\\ = 
B± = 0. For the 'thin boundary conditions', by means of the plasma shell, we have 
a regularization. The regularization parameter is the plasma frequency Q ^ oo. 
Hence, it would be useful to have a similar regularization for the 'thick boundary 
conditions'. One could imagine a plasma shell of finite thickness (which is a well 
investigated topic), allowing for a displacement ^ (in the bulk) not only parallel 
but also normal to the surface. This displacement field could be integrated out 
delivering a photon propagator which could be investigated in two different limiting 
procedures. The one would be to remove the regularization, i.e., to let Q ^ oo, 
which should deliver the 'thick boundary conditions'. The other way should be first 
to shrink the thickness of the plasma sheet to zero and after that to let — > cx). 
This should deliver the 'thin boundary conditions'. 
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Appendix 



In this appendix we consider a spherical plasma shell in order to give another il- 
lustration of the statement made at the end of section 4 that the excitations of the 
electromagnetic field and of the plasma have the same spectrum. We start from 
Eq. (125!) which is the kernel of the action of the displacement field after integrating 
out the electromagnetic field and which be means of Eq. fl26|) defines the equation of 
motion of the displacement field. First of all we have to consider the current j^{x), 
(jl]). The gradient entering jo contains derivatives in the tangential directions only, 
hence with Vj = riidr — -{n x L)i where n is the normal to the sphere and L is 
the orbital momentum operator, we get jo = e6{r — R)^{n x Further it is 

meaningful to introduce polarizations for the displacement field, 

^/'^ = L.^, = {n X L).^, (101) 



which form a basis. The corresponding expansion of the displacement fields reads 

2 
s=l 

We will see that the polarization s = 1 couples to the TE-mode of the electromag- 
netic field and s = 2 to the TM-mode. Now we consider the coupling with the 
electromagnetic field given in general form by Eq.([5]). In the spherical case we have 
instead of dHD now 

H|,^ = eSir - R) (^S^o^in x L)i - S^^do^ . (103) 
Being inserted into ([5]), the if^j appear in the combinations 

H^p!"^ = -edoE^^K (104) 
where we introduced the notations 

4" = £f = ( E r. ) ^ (105) 



fj. fj, 

with F = ^/^df+U/R^. These vectors E^'^ are orthogonal and normalized like 
the corresponding ones in (l62l) in the flat case. We note that these vectors up to 
differences in the notations coincide with those introduced in [T3]. 

The kernel ^K^, defined in Eq.f l25l) . carries indices like the displacement vectors 



^i. We project it onto the polarizations fllOll) and define 

(5^0)W^pt(^) (S^O)^^ pis) ^ (^Qg) 

Here we took into account diagonalness which holds for obvious symmetry reasons. 
Inserting into ( l25l) we get 

(5^0) (^) ^ ^0 _ (ij^^pWy D,^{z,z') {H^.Pf) . (107) 
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Here we took into account that = — mc^Q, Eq. fllUI) . remains unchanged and that 
is the free space photon propagator (122|) which we represent now in the form 



with 



I ^ ,./=o(.o-4) ^ Yi^iQ) d,{r, r') Y,l{n') (108) 



di{r,r') = ikoJi{hr^)hf\k^r^) . (109) 

We note that we use Lorentz gauge (a = 1) but any other choice would result in the 
same expressions for {^K^^ Traced back to (jl]) this is due to current conservation. 

Now we need to simplify fll07p . We define the corresponding quantities in mo- 
mentum representation, 

(5^0) (^) ^jdh ^..o(.o-4)^ YU^) {s^^f^ Yr^m (110) 



(note that [^K^Y"^^ is defined on the sphere). For the polarization s = 1 this i 



IS 



simple because the momentum operator L commutes with the propagator and we 
get simply 

{'K')12 = mkl + e'RXdi{R,R) 

^ mklgf\ko) (111) 

with 

g\'\ko) = {l + nR'diiR^R)) ■ (112) 

The corresponding calculation for the polarization s = 2 is a bit more involved 
because the operator (n x L) does not commute with the propagator. In the basis 
of the orbital momentum eigenfunctions we have to consider 

< lm\E^f g^uD{x - x')Ef'^\lm > . (113) 

The corresponding calculation was carried out in [T3j using known formulas from 
quantum mechanics and the result is 

(■^jsfO^j^) = mkl + e^^ (^^di{R,R)- <lm\{nx L)l D{x~x'){nx L)i\lm>^ 
^ mklgf\ko) (114) 

with 

9i [ko) = 1 " [r^^^ ' ' ~ ° 27Tl ) ■ ^ ' 

Using the recursion relations for the Bessel functions the last line can be rewritten 
in terms of derivatives of the Riccati-Bessel functions, 

gf'ik^) = 1 + ^ jl{koR)h[{koR) . (116) 
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We finish this part with the remark that the functions Qi (ko) are the same as the 
Jost functions for the corresponding scattering problem for the electromagnetic field 
in a spherical plasma shell. The corresponding scattering phase shifts are given in 

Returning to Eq. fl25|) it is clear that after expanding the displacement fields in 
the orbital momentum basis the equations of motion become algebraic ones like 
in the case of fiat plasma sheet and that the spectrum is defined by the zeros of 
the functions gl'^\ko)- In this way, for a spherical shell, we demonstrated that the 
excitation spectrum of the electromagnetic field and of the displacement field are 
the same. Further, since it is known that these functions do not have zeros for real 
ko we come to the conclusion that the displacement field does not have plasmon 
excitations. We conclude this appendix with the remark that the relation between 
the spectra of the excitations of the electromagnetic field and the displacement field 
is, of course, more involved. What we showed here is merely that the excitations of 
the displacement field are also excitations of the electromagnetic field since the zeros 
of ^k^ fl25l) are poles of fl28l) . But the electromagnetic field has more solutions, 
namely the scattering states, which do not have a correspondence in the spectrum 
of the displacement field. 

We remark, that of course, since we were looking for solutions with real u, we 
have in mind excitations not decaying in time. Further we would like to mention 
that the absence of excitations found here is not in contradiction to the excitations 
found in [3] since those were obtained in an approximation which was nonrelativistic 
from the outset. 
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